Increasing entanglement through engineered disorder in the random Ising chain 
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The ground state entanglement entropy between block of sites in the random Ising chain is studied 
by means of the Von Neumann entropy. We show that in presence of strong correlations between the 
disordered couplings and local magnetic fields the entanglement increases and becomes larger than 
in the ordered case. The different behavior with respect to the uncorrelated disordered model is 
due to the drastic change of the ground state properties. The same result holds also for the random 
3— state quantum Potts model. 
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Entanglement plays a central role in modern quan- 
tum mechanics. Having been regarded for long time 
as the root of the incompleteness of quantum mechan- 
ics and the source of several paradoxes [H, with the ad- 
vent of Quantum Information (QI) theory it has received 
renewed attention, and acquired the status of a funda- 
mental resource for QI processing 0]. In the context of 
strongly correlated quantum systems entanglement arises 
in a natural way [3|, and can in fact be regarded as 
a conceptual bridge between condensed matter physics 
and QI theory. In particular, much attention has been 
recently devoted to the ground state entanglement be- 
tween two blocks of spins in one dimensional spin chains 
i, H, i, 0, S S 0, [H, [H, [B, P . By means of both 
analytical proofs and numerical simulations, this quan- 
tity has been shown to be logarithmically divergent with 
the length of the block when the spin chain is critical. Re- 
markably, the prcfactor of the logarithm is proportional 
to the central charge of the corresponding conformal field 
theory (CFT) associated with the spin model at hands 

The average block entropy has been also studied in the 
context of random critical spin chains [H, [13, [H, [H, HO, 
HHj III this case, on the one hand works on random 
Ising and XXZ spin 1/2 chains [l^, suggest that 

the block entropy still grows logarithmically with the size 
of the block but with a smaller prcfactor of the logarithm 
with respect to the ordered model, being "rcnormalized" 
by a factor In 2. On the other hand, the block entropy in 
random quantum Potts chains with spin dimension d > 2 
has been studied in [l^, where a strong-disorder renor- 
malization group (RG) analysis shows that the average 
entropy still diverges logarithmically but that the pref- 
actor of the entropy is not proportional to the one of the 
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corresponding pure model and becomes larger than the 
latter when d > 41. This provides evidence that disorder 
can increase entanglement with respect to pure models, 
contrary to naive expectations. The same analysis have 
been performed in aperiodic spin 1/2 chains and again 
it has been shown that the prefactor of the entropy log- 
arithmic scaling is not simply the "renormalized" factor 
of the homogeneous model [23| . Therefore at present the 
possibility of associating an effective central charge to a 
non-homogeneous system is unclear. 

In this work we study a random critical Ising spin chain 
j2l | where the (nearest neighbor) couplings and the lo- 
cal transverse magnetic field are drawn from the same 
probability distribution and in addition share a certain 
(tunable) degree of correlation. We compute the average 
block entropy for this model and show that the prefac- 
tor of the entropy in a class of random correlated chain 
is larger than in the homogeneous case. An analogous 
result is found to hold also for larger dimension as we 
demonstrate in the random correlated quantum Potts 
model with d = 3. 

We consider a spin chain with open boundary condi- 
tions and with d states |0), |1), . . . , jd — 1) per lattice site, 
described by the Hamiltonian (2^ 
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Here L is the length of the chain, P represents the ladder 
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operator r|s) = |(s-Hl) mod d), (/c| S"^ 
{k,k' — 0, . . . , d — 1) and is the hermitian conjugate 
of 5". 

For d = 2, 5^ = CT^ (with a'-^'-y''^ the Pauh matrices) 
and the model reduces to the (random) transverse field 
Ising model 



L-l 



H2 = - E J^'^i'^i+i 



L 



(2) 



In the pure case Ji ~ hi ~ I while in the random case Ji 
and hi are positive random numbers drawn from a joint 



distribution Va{J, h) defined in the interval [0, 1]^^ such 
that the marginal probability distributions 



P(J) 



dhVa{J,h) 



P{h) = / dJVa{J,h) 
Jo 



(3) 
(4) 



are both uniform distributions in the interval [0, 1]. This 
ensures that the random models we consider are criti- 
cal [111- The correlation coefficient characterizing the 
joint distribution Va is defined as a = ((J — {J)){h — 
i^))) /i'^J'^h) where ctj (ct/j) is the standard deviation 
of J (h). In the following we will consider the interval 
< a < 1 whose endpoints correspond to independent 
and perfectly correlated couplings and transverse fields, 
respectively. 

We study the ground state entanglement between a 
block of the first i spins and the rest of the chain. This 
is quantified by the von Neumann entropy of the reduced 
density matrix of the block: 



Se = -Trpf log2 /Of 



(5) 
(6) 



with \iPg) the ground state of the Hamiltonian. In the 
random case we average over N independent disorder re- 
alizations. 

For homogeneous models, CFT calculations show that 
the entropy scales as Q: 
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where k is equal to the central charge c of the correspond- 
ing theory and A is a non-universal constant . The 
central charge corresponding to the homogenous quan- 
tum Potts model depends on the local site dimension d 
as 



^d-1 
'd+2' 



(8) 



[l9j. For general non- homogeneous models, as said be- 
fore, the entropy still grows log arithmically but with a 
different prefactor k^c [13, [11 HI • 

The d = 2 model can be solved analytically in both the 
homogeneous and disordered case as it reduces, as no- 
ticed earlier, to the transverse field Ising model [H, [13]: 
one first performs a Jordan- Wigner transformation that 
maps the spins into fermions, then using a Bogoliubov 
transformation, whose coefficients are found via an ex- 
act diagonalization, the system is mapped into non- 
interacting fermions. Finally, the eigenvalues of the re- 
duced density matrix and thus the entanglement can 
be determined from the two-point regular and anomalous 
correlation functions as shown in [28| . 

For systems with higher dimensionality one cannot re- 
sort on any analytical tools, and in order to study the 
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FIG. 1: (color online) Block entropy (in bits) averaged over 
N = 5 X 10^ configurations for an Ising spin chain of L = 200 
sites in the homogeneous (squares) and perfectly correlated 
random (circles) cases. Best-fit lines with expression ((Tjl are 
also shown. 



3— states quantum Potts model we turn to the finite size 
density matrix renormalization group (DMRG) method 
t3Q., 2L, 2M- The strategy behind the DMRG algorithm 
is to construct a system block and then recursively en- 
large it, until the desired system size is reached. At ev- 
ery step the basis of the corresponding Hamiltonian is 
truncated, so that the size of the Hilbert space is kept 
manageable as the physical system grows. The trunca- 
tion of the Hilbert space is performed by retaining the 
eigenstates corresponding to the m highest eigenvalues 
of the block reduced density matrix. Several sweeps of 
the finite system DMRG are performed to increase the 
accuracy for non-homogeneous chains. 

d=2 results: Ising model ~ In FigU] we show the re- 
sults for the block entropy evaluated for the homogenous 
and perfectly correlated random (a = 1, = 5 x 10'^ 
configurations) Ising chains of total length L = 200. Fit- 
ting the numerical data with expression ([7]) we obtain 
for the homogeneous system = 0-50 ± 0.01, in perfect 
agreement with Eq.([8]) which gives C2 = 0.5 for the homo- 
geneous chain with d = 2. We obtain ^2^^ = 0.61 ± 0.01 
for the perfectly correlated a = 1 case: Notice that, dif- 
ferently from previous studies on random spin one half 



chains, ^2^ is larger than the prefactor for the homo- 
geneous case resulting (in the thermodynamic limit) in 
an increment of the entanglement in the disordered chain 
with respect to the homogeneous case [33| . 

An important question that ought to be studied is the 
stability of this behavior with the correlation a shared 
by the couplings Ji and the transverse fields hi. In Fig[2] 
we show the prefactor of the average entropy for dif- 
ferent chain lengths (L = 50, 70, 100) as a function of 
the correlation a. The results show that for small a we 
confirm the analytic result of [iTj apart from finite size 
errors. Indeed, a finite size scaling (see inset of Figl^]) 
gives ^2^° = 0.30 ± 0.04 quite in agreement with the 



3 




0.2 0.4 0.6 0.8 

a 



FIG. 2: (color online) Prefactor of the entropy versus the 
correlation coefficient a of the J and h couplings for different 
chain lengths (L = 50, black circles; L — 70, red squares; L = 
100, blue diamonds) and N = 10'* configurations. The solid 
and dashed line gives Kg = C2 In 2 ~ 0.35 and Kg = C2 = 0.5 
respectively. The inset shows the finite-size scaling of the 
prefactor with best-fit lines (solid, a = 0; dotted, a — 0.95; 
dashed-dotted, q = 1). Extrapolation to the thermodynamic 
limit 1/L ^ gives kJ=° = 0.30±0.04 and kJ=i = 0.60±0.01. 



expected result ^2^° — C2ln2 « 0.346; in the perfect 
correlated case one gets instead ^2^^ = 0.60 ± 0.01 in 
the thermodynamic limit. For intermediate values of a 
our results show the departure from the completely cor- 
related random Ising model and for strong enough cor- 
relations (a > a* ^ 0.9) the prefactor become larger 
than the prefactor for the pure model = 0.5. This re- 
sult is a clear evidence that it is possible to increase the 
entanglement adding disorder to an homogenous system 
already for small local dimension d = 2. Notice that the 
scaling of K2 with the system size is faster for stronger 
correlations. Indeed, the perfectly correlated case, a = 1, 
seems at convergence already with L = 50 while the un- 
correlated case is still far for the thermodynamic limit at 
L = 100. This clearly reflects the largest space of the 
possible configurations of the uncorrelated case. How- 
ever, the finite size scaling suggests that in the whole 
range a* < a < 1 we have > '*2 ~ '^2 in the thermod- 
inamical limit. 

A deeper insight on the ways the a > a* case differs 
form the a ~ one, can be gained by determining the 
probability distribution P of the block entropy for a fixed 
block length in the two cases. This procedure will in fact 
unravel the presence of long distance effective quantum 
correlations in the chain ground state. The resulting dis- 
tributions for L ~ 100, ^ = 50 and = 5 x lO'^ are shown 
in Fig. [3l From these results it emerges that the ground 
states for the two models are structurally different: for 
a = a single narrow peak at 5*50 = 1 denotes the signa- 
ture of the random singlet-like phase [13, [11] ; for a = 1 
the peak at S50 = 1 is wider indicating the occurrence 
of a more complex structure. This is a clear signature 
that the clustering procedure described in [13] breaks 




FIG. 3: (color online) Probability distribution for the von 
Neumann block entropy for a block of length ^ = 50 of an 
L — 100 chain for N = 5x 10^ configurations. Left panel: a = 
(random uncorrelated case); right panel: a = 1 (random 
perfectly correlated case.) 



down. Indeed, the strong disorder renormalization pro- 
cedure rely on the renormalization of the strongest bonds 
or local field, obtaining either a ferromagnetic cluster or 
a frozen spin. However, if the magnetic field and the cou- 
pling are always of the same order this renormalization 
procedure may not be applied. It seems there exists a 
crossover between the two extreme regimes, namely, the 
random-singlet phase and the delocalized one (see Fig. 
IS). 

d=3 results: Potts model - In Fig[3]we show the block 
entropy for L = 70 in the homogeneous and perfectly cor- 
related random cases for d = 3, evaluated through a finite 
size DMRG algorithm using m = 30 and = 3 x 10^ 
configurations. With these parameters the maximum er- 
ror in the entropy of a single configuration is less than 
5 X 10~^. The fit of the data for the pure model according 
to expression ([7]) gives k| = 0.83 ± 0.02, in good agree- 
ment with value predicted by the CFT formula of Eq. ([5]) 
C3 = 4/5 = 0.8. In the disordered case we find a similar 
picture as for the d = 2 case: in the perfectly correlated 
random case we obtain a bigger value of the prefactor 
tq=^ = 0.91 ± 0.02 > k\. For a = it is lower than 
that for the homogeneous model Kg^^ = In 3/2 w 0.55 
(computed in [T9| via the strong disorder renormaliza- 
tion group). 

Conclusions - In this paper we have considered a class 
of random critical spin chains described by the Hamil- 
tonian of Eq.([T]). We concentrated on two models: the 
random Ising model with d = 2 and the quantum Potts 
model with d = 2>: the couplings J and the local trans- 
verse magnetic fields h are random but share a given 
amount of correlations. Our analysis shows that when- 
ever J and h are correlated the prefactor of the Von Neu- 
mann entropy is always greater than the prefactor for the 
uncorrelated case. Moreover, when the correlations are 
above a certain threshold a* the prefactor of the ran- 
dom correlated model is larger than the prefactor for the 
homogenous chain. It would be interesting to confirm 
our numerical results applying a strong disorder renor- 
malization procedure. In this work we first demonstrate 
how to increase entanglement by putting static disorder 
in systems with small local dimension. It is not needed to 
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FIG. 4: (color online) Block entropy (in bits) averaged over 
= 3 X 10^ configurations for a d = 3 Potts model of L = 70 
sites in the homogeneous (squares) and perfectly correlated 
random (circles) cases. Best-fit lines with expression (0 are 
also shown. 

consider systems with high local dimensions d > 41 as in 
[l^ : our results indeed shows that exists a class of disor- 
dered spin one-half models where the ground state entan- 
glement is larger than in the correspondent translational- 



invariant ones. Notice that this random correlated model 
could be realized experimentally in engineered quantum 
systems, e.g. in optical lattices, ion traps or arrays of 
Josephson junctions. We conclude observing that the 
correlated random models studied in this work lay out- 
side the random singlet-like phase for which a generalized 
c-theorem holds. It is thus ambiguous the association of a 
renormalized central charge to general non-homogeneous 
systems. 
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